


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1964 


Digital computer solution of linear and 
nonlinear ordinary differential equations 


Cooper, David S.; Knepper, Donald E. 


Monterey, California: U.S. Naval Postgraduate School 
http://ndl.handle.net/10945/12012 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 


1964 
COOPER, D. 





DIGITAL COMPUTER SOLUTION OF LINEAR AND 
NONLINEAR ORDINARY DIFFERENTIAL EQUATIONS 


DAVID S. COOPER 
DONALD £, KNEPPER 


apes LIBRARY 
ALAN STGORADUATE SCHOO 
MONTEREY CA 93943-5101 i. 


LIBRARY 
U.S. NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 




















DIGITAL COMPUTER SOLUTION 
OF LINEAR AND NONLINEAR 


ORDINARY DIFFERENTIAL EQUATIONS 


et ce 


David S. Cooper 
and 


Donald E. Knepper 





DIGITAL COMPUTER SOLUTION 
OF LINEAR AND NONLINEAR 


ORDINARY DIFFERENTIAL EQUATIONS 


by 
David S. Cooper 
Lieutenant, United States Navy 
and 
Donald E. Knepper 


Lieutenant, United States Navy 


Submitted in partial fulfillment of 
the requirements for the degree of 


MASTER OF SCIENCE 
IN 
ELECTRICAL ENGINEERING 


United States Naval Postgraduate School 
Monterey, Californis 


1964 








= or and waty e : 
ry at tf @ “Reg ytys 7 4.50) A bran 
DUDLEY XNOX LIBRARY 
, f A 

fA 


LPostce 


MONTEREY Cc 


DIGITAL COMPUTER SOLUTION 
OF LINEAR AND NONLINEAR 
ORDINARY DIFFERENTIAL EQUATIONS 
by 
David S. Coopér 
and 


Donald E. Knepper 


This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 
ELECTRICAL ENGINEERING 
from the 


United States Naval Postgraduate School 








ABSTRACT 

A digital computer program which will provide accurate and stable 
solutions of Linear and nonlinear ordinary differential equations is 
described, as are subprograms which simulate many of the nonlinearities 
occurring in feedback control systems, The program employs fourth 
order Runge-Kutta numerical integration with automatic error checking 
and interval modification. Specialized coding sheets assist in the 
preparation of input data, and built-in print and graph output routines 
provide a permanent record of input equations, parameter values and 
output data for each solution. The program is flexible, yet it can 
be used by a person with only a rudimentary knowledge of FORTRAN 


programming. 
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Le Introduction. 

The analysis of scientific problems frequently leads to a set of 
simultaneous ordinary differential equations which describe the phenomenon 
being studied. Analytic solution of these equations is often impracti- 
cable and analog or digital simulation becomes necessary. A digital com- 
puter program which will solve ordinary differential equations has been 
written in FORTRAN 60 for the CDC 1604 Digital Computer. It includes 
FUNCTION subprograms which cover many of the nonlinearities occurring 
in feedback control systems | 7] . The program has been named SIDES, 
an acronym for SImultaneous Differential Equation Solver. 

Digital simulation should be considered as a valuable supplement 
for the analog computer, but not as a replacement. Selection between 
analog and digital simulation methods should be based on an evaluation 
of the problem to be solved: the complexity of the equations, the magni- 
tude of the included parameters and the required accuracies of the solu- 
tions. Factors which favor the selection of digital simulation, and 
more particularly, SIDES, are: 

a. The necessity for number scaling is eliminated. 

b. Solution accuracies are achieved which are unattainable 
with an analog computer. 

c. Printed and graphical outputs provide a permanent written 
record of input equations, initial parameter values and output data for 
each solution. 

d. Nonlinear operations such as squaring, forming trigono- 


metric and exponential functions, and multiplication are easily performed. 








e. Special functions such as hysteresis and time delay are 
available or can be readily generated as ordinary FORTRAN SUBROUTINE 
or FUNCTION subprograms. 

Some of the limitations are: 

a. Inherent errors are due to the sequential operation of the 
digital computer in contrast to the simultaneous processing of signals 
performed by the analog computer. 

b. The effect of parameter variations on problem results cannot 
be immediately observed. (This limitation could be partially eliminated 
by the incorporation of a data display unit with the CDC 1604). 

Existing digital simulation programs vary considerably with regard 
to the input language and the integration techniques used. Stein, Rose 
and Parker[5| have devised a system which accepts as input the descrip- 
tion of a "patching" diagram prepared for the Electronic Associates PACE 
Analog Computer. FORTRAN is used as an intermediate in the compilation. 
The DYSAC program written by Hurley [2] is in symbolic machine language 
(CODAP). It contains a supply of digitally simulated analog computer 
components which are interconnected by the "patching" information sup- 
plied by the user. Integration is performed by the fourth order Runge- 
Kutta method with a fixed integration interval. Nordsieck [4] has 
developed a method for automatic digital computer integration which is 
equivalent to a reformulation of the method of Adams. 

The SIDES program to be described here, was designed with particular 
emphasis on user convenience and solution accuracy. Input and output 
formats are built-in and integration is performed by the fourth order 
Runge-Kutta method with automatic error checking and interval modification 
as developed by Anderson [1] . A listing of the SUBROUTINE and the FUNCTION 


2 


subprograms which have been developed for use in conjunction with SIDES 


is given in Appendix B. 





om The SIDES Input Routine. 

Calling program. SUBROUTINE SIDES will normally be called from the 
FORTRAN Library tape by a calling program which must be of a prescribed 
form. To simplify the writing of this program, a SIDES Coding Sheet has 
been prepared (Appendix A) which reduces most of the job of programming 
to the simple task of filling in blanks and lining out undesired options. 


The calling program deck must be of the following form: 


.. JOB* JONES J BOX 113 (Job ident.) 
PROGRAM PROBLEM (Arbitrary name) 
DIMENSION X(50), XDOT(30), C(50) (Standard card) 
G50) “="r, 0 , " an 


1 CALL SIDES (T,X, XDOT,C) - ie 
(Equations written 
in FORTRAN language) 
GO TO 1 (Standard card) 
END i 
END i a 
This deck is followed by a blank card and the data and option cards. 
A typical program, together with its output, is included in 
Appendix A, 
Equation statements. While equations for different problems vary, 


they can be manipulated into the standard first order form: 


dx, : 
f= £ (xy o Moe ees X 5 t) i aoe) 


dt 


in which the independent variable must be designated by t and the 





dependent variables by Xo The number of such equations is designated 
by n with nS 30. Since only the x and t variables can be output, addi- 
tional Xo 0 < ix 48, have been made available. If some quantity other 
than one of the variables is desired as output, it must be equated to 
one of these x,. 
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A set of constants, c l= 49, has also been made available so 


i? 
that parameters may be introduced or varied by data cards. 

Any legal FORTRAN statement or technique may be used in writing 
these equations, provided that a transfer of control does not prevent the 
processing of all relevant statements during each step of numerical 
integration. All of the FORTRAN Library functions are available for 
appropriate inclusion in the equation statements. Some additional 
FUNCTION subprograms which have been developed for use by SIDES are 
described below. These and other subprograms may be appended to the 
calling program in the conventional manner. They may also be added to 
the library tape if appropriate. 

Data and option cards. The seven sections of input data required 
to specify a problem for SIDES are enumerated in Table I. All sections 
are straight-forward and their purpose fairly obvious. 

If several solutions of the same equation are required, it is 
necessary to specify the number of runs and to supply a data deck for 
each run. For each run after the first, only the cards of sections 2 
through 7 (Table I) are required. Information can be retained between 
runs by choosing the appropriate HELD option. If it is desired to sub- 
divide a long solution into two or more successive runs, the final condi- 


tions of one run may be used as initial conditions for the next by making 


use of the FCS HELD command, 
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TABLE [I 


DATA REQUIRED FOR A PROBLEM 


DATA SECTION 


Problem description 


Coefficients 


Initial conditions 


Time 


Integration step 
size 


Data print-out 


Graph output 


CONTROL OPTIONS 


1— 30 first order differential 


equations 
1—9 runs 


Values zeroed 
Values set in 
Values held 


Values zeroed 

Values set in 

Initial or final condition 
values held 


Initial and final time 
given 
Values held 


Step size given 

Step size computed and 
automatically varied 
to maintain a pre- 
scribed precision 

Step size held 


1— 8 variables printed 
1—> 99 steps of integration 
between print-outs 

Data held 
No print-out 


1—>5 single-curve graphs 

1—»>5 curves on a single 
graph 

1—> 99 steps of integration 
between graph points 

Data held 

No graph output 





Oe The SIDES Output Routine. 

All job identification, together with the run number and a record 
of all input data, is automatically output for each run. A maximum of 
eight quantities may be printed-out during integration. An automatical- 
ly scaled graph output is also available, providing a maximum of either 
five single-curve graphs or five curves on a single graph. 

A maximum of 10,000 integration increments, 300 lines of print-out, or 
900 graph points (per curve) can be generated during any one run. The run 


will be terminated if any one of these conditions is exceeded. 
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4, Error Indications. 

The SIDES program will be terminated when all runs called for are 
completed or if an error is detected in the calling program or in the 
data and option cards. SUBROUTINE SIDES contains 38 error stops, each 
preceded by a diagnostic print-out. These error stops refer to diffi- 
culties such as incorrectly punched cards, missing cards and improper 
option selection. All error stops and print-outs inherent to the basic 


FORTRAN Compiler also apply. 





Di, The Integration Method. 

The operation of SIDES is centered around the fourth order Runge- 
Kutta method of numerical integration, The choice of integration algo- 
rithm is necessarily a compromise and Runge-Kutta was chosen primarily 
because of the ease both in starting solutions and in changing the 
integration increment during solutions. The use of the method developed 
by Anderson [1] to compute the truncation error has eliminated one of 
the main shortcomings of Runge-Kutta as used in other programs [6| . 

On this basis, the integration interval can be automatically varied to 
ensure stability and to keep the error during each step within prescribed 
bounds, without excessive lengthening of the computing time. The allow- 
able truncation error per unit time for each of the dependent variables 
was chosen to be: 


ne 10 * (variable magnitude ) 


i total time 
The variation of the step size is bounded by arbitrary limits of 10°? 
and 10° "times the total time. 
In the event that the user wishes to prescribe a fixed step size, 


the Runge-Kutta integration can be processed in either one or two 


segments, each with a different step size. 








Gr Accuracy Considerations. 


Since the accuracy of the Runge-Kutta method of numerical integra- 


tion is a function of the size of the integration increment, use of a 


digital computer allows solution accuracies that in most cases far ex- 


ceed existing engineering demands. 


The following two examples are 


given to indicate the effect of integration interval on solution accuracy. 


The interval was held constant for each run. 


Example 1! 


Differential Equation 


Explicit Solution 


Integration Interval 
0.03 


Example 9? 


Differential Equation 


Explicit Solution 


Integration Interval 
0.005 


vid ees 


3.8x107" 


3.8x107- 


3.8x107- 


1 _ 2 
1+x 


y = (14x)? 


y at x5e09 


9,4%10° (at 


eee 2.59) 


8.1x10> 


8.1x10° 


(yo0= 1, xa) 


Error at x = 5 


x10° } 


3x10"? 


iealtsy 


(y .= 1 Ko= 0) 


Error at x = 5 


sb ol (at x=2.9) 
Done 


4x10° 


Sider W. E. Numerical solution of differential equations. 
John Wiley & Sons., 1953: pp. 20-21. 


2 tbid., PpPpe 74-75. 
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Although Example 2 is intentionally unfavorable for solution by 
Runge-Kutta integration, it is obvious that high accuracies can be 
achieved, although at the expense of additional computing time. 

One of the main advantages of the computed step size integration 
method is that solution accuracy and computing time are not seriously 
impaired by problem discontinuities, This is an important factor to 
the engineer who must analyze systems containing switching devices or 
step-wise parameter fluctuations. Solution accuracy in SIDES is deter- 
mined by the estimated magnitude (power of ten) of the dependent vari- 
ables, X(j), 1< j<n. An examination of Example 3 indicates that 
this estimated magnitude can often vary widely before solution ac- 


curacy suffers. 


Example 3 





where R is of the form: 
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VARIABLE MAGNITUDE 


SETTING 





ERROR AT 


T=2 


dsc 


5x10 +? 


5x10"? 


1x10’ 


2x107° 


7x10" ° 


7x10°° 


7x10" ° 
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ERROR AT 


T= 4 





ie Special FUNCTION Subprograms. 

One of the stated advantages of digital solution was the relative 
ease with which many nonlinearities can be simulated. The develop- 
ment of subprograms to simulate functions such as hysteresis and time 
delay, however, presented special problems in relation to their use with 
Runge-Kutta integration with variable step size. Since FUNCTION simula- 
tion of these nonlinearities requires past information for its proper 
operation, it was necessary to modify the basic SIDES program so that 
there could be a free interplay of information between it and the 
particular FUNCTION subprogram. Utilizing the fact that the COMMON 
statement, as applied to FORTRAN 60, will cause values to be entered 
in known cells, it was possible to establish this needed correspondence 
between the two separately compiled subprograms. 

FUNCTION subprograms have been written to simulate time delay and 
two and three position relays with hysteresis. These subprograms are 
available for use in conjunction with SUBROUTINE SIDES. A brief descrip- 
tion of each FUNCTION is included in Appendix A anda listing of each 


FUNCTION subprogram appears in Appendix B. 


is. 








8. Recommendations and Conclusions. 

SUBROUTINE SIDES has proven to be an extremely easy program to use. 
The SIDES Coding Sheet enables the user to prepare his program with a 
minimum of reference to additional explanatory material, Due to the 
complexity of the input program, it is believed that the coding sheet 
is almost a necessity. Continued use of this sheet will probably point 
out needed changes to increase its effectiveness. It is already antici- 
pated that a reduction in the data input for successive runs is both 
feasible and desirable. 

Although fourth order Runge-Kutta numerical integration with auto- 
matic error checking and interval modification is accurate, stable and 
relatively rapid, it is recommended that other integration methods in- 
cluding predictor-corrector types, and third and fifth order Runge- 
Kutta be evaluated, 

The FUNCTION subprograms developed for use with SIDES consist of a 
few of the nonlinearities encountered in the analysis of feedback con- 
trol systems. Consideration has been, and should be further given to 
the development of differentiators, noise generators and function genera- 
tors. 

SUBROUTINE SIDES is intended to provide a simple but moderately 
flexible means for the solution of differential equations. The inte- 
gration technique provides a high degree of accuracy and stability and 
the prepared coding forms simplify the programmer's task. [It is be- 
lieved, therefor, that SUBROUTINE SIDES is a valuable and essential 


supplement to the analog computer. 
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SUBROUTINE SIDES PROGRAMMING INSTRUCTIONS 


SUBROUTINE SIDES is intended to facilitate the solution of a class 
of problems involving linear and nonlinear differential equations. It 
was written in FORTRAN 60 for the CDC 1604 Digital Computer. The sub- 
routine will normally be called from the FORTRAN Library tape by a 
program which is being compiled and executed under monitor control. 
This calling program and the associated seven data sections required 
by SIDES for the solution of the problem are supplied on 80 column 
punched cards. Utilization of the SIDES Coding Sheets will greatly 
simplify the programmer's task and serves as a check that problem data 
and control commands are properly positioned and in the correct order. 
A description of the functional sections of this form and instructions 
for their use are given below. In what follows, reference should be 
made to the copy of the coding sheets at the end of this appendix. 
These have been filled out for a typical problem and the resulting 


computer output is also appended. 


1. The Calling Program, 

This section contains the problem equation statements which will 
probably cause the inexperienced programmer the most difficulty. 

The first card contains the job identity, normally just the pro- 
grammer's name and his computer facility box number. 

The second card contains the program name. Any combination of 
seven letters or less is acceptable. 

The next three cards are standard and must be punched as shown on 
the coding sheet. 
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Following these cards is the statement of the differential equa- 
tions to be solved. These equations must be written in the form of n 


simultaneous first order equations, where n must not exceed 30. Thus, 


XDOT(1) = £,(XC1), X(2),...,X(n),T) 


e e e e e e @ e 


XDOT(n) = £ (XC), X(2),..+5X(n),T) 
where, 

ah = the independent variable 

X(j) = the jth dependent variable 


XDOT(j) = the first derivative of X(j) with 
respect to T 

X(j) may be identified as the output of the jth integrator in the 
corresponding analog computer setup, while XDOT(j) may be identified as 
the input to that integrator. Notice, however, that the sign reversals 
inherent in the analog operations do not occur here. 

It is often desirable to define other functions of the above vari- 
ables, either so that the functions, oe can be built up in more than 
one step, or so that a quantity other than one of the variables can be 
printed or graphed. Although any acceptable FORTRAN variable name may 
be used to define these functions, they must be equated to an X(j), 
n<j< 48, if they are desired as output. 

Finally, if several solutions are contemplated with changes in some 
parameters, the constants C(1) through C(49) may be introduced. Their 
values may then be entered later on data cards. The constant C(50) 


must never be used except as in the calling program, above. 
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Block diagram manipulations are often useful in the formulation 
of the desired first order differential equations. The following 
example shows how a simple pole-zero transfer function may be reformed 


and expressed in FORTRAN equation statements. 


- s+6 





s+3 


3 
pee”: S+6 s+6 


R = C(1)4T 

ERROR = R - X(2) 

X(2) = ERROR - X(1) (permits print-out or graph 
XDOT(1) = 3.0*(ERROR - 2.0*X(1)) a 


X(3) = ERROR (permits print-out or graph 
of the ERROR) 


Any legal FORTRAN statement or technique can be used in the call- 
ing program, provided that a transfer of control does not prevent the 
processing of all relevant statements during each step of numerical 
integration. Improper transfer of control has proven to be one of the 
most prevalent types of programming errors made by users of SIDES. In 
the variable step size mode of integration, the equation statements are 


frequently evaluated with too large a step size and consequently too 
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large a truncation error. These values are then automatically reject- 
ed by the program and the equations are re-evaluated with a smaller 
step size. This process is repeated until the truncation error is less 
than the prescribed amount. Care must therefor be taken to ensure that 
values which have been set after a transfer of control on an unaccept- 
able pass are reset correctly on the accepted pass. 

The COMMON statement must not be used unless the programmer is 
completely familiar with the SIDES program. All of the FORTRAN Library 
functions are acaba for inclusion in the equation statements. Addi- 
tional FUNCTION subprograms applicable to common control system non- 
linearities have been developed for use in conjunction with SIDES 
and will be discussed later. If the user desires to write his own 
FUNCTIONS or SUBROUTINES, they may be added to the calling program in 
the conventional manner. Comments are permissible throughout the call- 
ing program if a C is placed in column one. 

The calling program must be terminated by: 

GO TO l 
END 
END 


blank card 


2a The Data and Option Cards. 

C Values. One of the options in parenthesis must be Lined out. 
The values of the coefficients, C(j), not simply zeroed or held, are 
placed, with decimal point, in any order, on the next three lines. 


Each coefficient is preceded by its subscript, j, right justified, in 
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the appropriate block. Unused lines are lined out and the total 
number of lines (Q0->3) containing values is noted on the option 
line for every run. 

X Values. All but one of the options must be lined out. On 
any run after the first, either the initial or final conditions of 
the preceding run may be held and automatically set as the initial 
conditions of the new run. The values are entered in the same manner 
as the C values, lining out the unused Lines. 

Time. On the initial run, the initial and final problem times 
are entered, with decimal point, and the HELD command is lined out. 
On subsequent runs either option may be chosen. The final time has 
an upper limit of 10,000 time units. 

Integration Step Size. All but one of the options must be lined 
out. 

If the step size(s) are given, their value(s) are entered, with 
decimal points, on the next line with the corresponding termination 
time(s). The last specified termination time must be the same as the 
final time value entered in the preceding section. Line out the next 
two lines, 

If the step size is to be computed, the line for step size data is 
lined out and the order of magnitude, N, (X(j) > io’, -9<N <99) of 
each dependent variable, X(j), 1<j<n, is entered, right justified, 
following the appropriate variable subscript on the next two lines. 
Unused variable subscripts should be lined out. If only one dependent 


variable appears in the problem, the second line is lined out. 
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If the HELD command is given, the following three lines are lined 
out. 

Data Print-out. All but one of the options must be lined out. 

If the GIVEN command is selected, the number of integration steps 
between print-outs must be entered in the appropriate block. Although 
a print-out at every 20th step will usually be sufficient, it should 
be kept in mind that the program will automatically stop at an upper 
limit of 300 lines of print-out. The next line must contain the sub- 
scripts, right justified, of the variables to be printed. For print- 
out, the independent variable, T, is specified by the subscript 50 and 
the integration step size by the subscript 49. On the next two lines 
enter the data print-out column headings in the appropriate block under 
the corresponding variable subscript. All three data cards must be 
present when this option is chosen. 

If the HELD or NONE commands are chosen, the following three lines 
are lined out. 

Graph Output. All but one of the options must be lined out. 

If graph output is needed, enter the desired number of graphs or 
curves in the appropriate block. On the next line enter the number of 
integration steps between plotting points. This should be kept small 
to ensure smoothness, remembering, however, that the program will auto- 
matically stop at an upper limit of 900 plotting points per curve. On 
the last line, enter in the designated blocks the variable subscripts 
of the X and Y ordinates, respectively, for graphs (curves) A, B, C, 

D, and E. Unused graph (curve) letters may be lined out. Since the 


first curve on the multi-curve graph sets the scale, the "largest" 
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curve must always be specified first. 
If the HELD or NONE commands are chosen, the last two lines are 


lined out. 


ais Special FUNCTION Subprograms. 

Certain FUNCTION subprograms have been developed for use in con- 
junction with the simulation of linear and nonlinear feedback control 
systems by SUBROUTINE SIDES. These subprograms are called into use by 
the normal FORTRAN techniques and may be used with either integration 
method of SIDES. The arguments may be expressed by variable names or 
numerical values, with decimal point. A brief description of each 


FUNCTION and its arguments follows. 


FUNCTION DELAY! (DELT, VAR, T) 

FUNCTION DELAY! will provide as output the designated variable'’s 
value, delayed in time by a prescribed amount. The output of the de- 
lay will be the initial value of the variable until such time as T 
exceeds DELT. Since this FUNCTION involves storage of past informa- 
tion, it can only be used once in a given program. For this reason, 
additional subprograms, FUNCTION DELAY2 and FUNCTION DELAY3 are also 
available. The arguments of each of these subprograms are: 

DELT The value of delay time. The length of delay must 
be greater than the integration step size and less 
than the sum of the previous 499 time increments or 
an error print-out will be given and the program 


will stop. 
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VAR The variable name of the quantity being delayed. 
T The independent variable, time. It must be represent- 


ed by T in the calling statement. 


FUNCTION DZONE (XINPUT, HALFD, GAIN) 


FUNCTION DZONE simulates dead zone or free play. 
OUT 





Its arguments are: 


XINPUT The variable name of the input quantity. 


HALFD The absolute value of one half of the dead zone 
or band. 
GAIN The value of the gain of the nonlinearity. 


FUNCTION REL2 (XINPUT, OUTPUT) 


FUNCTION REL2 simulates an ideal two position relay. 
OUT 


OUTPUT 


IN 
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Its arguments are; 
XINPUT The variable name of the input quantity. 


OUTPUT The absolute value of the output of the relay. 


FUNCTION REL2H1 (XINPUT, HALFD, OUTPUT) 


FUNCTION REL2H1 simulates a two position relay with hysteresis. 
Since it involves the storage of past information, it can only be used 
once in a given problem. For this reason, additional subprograms, 
FUNCTION REL2H2 and FUNCTION REL2H3 are also available. 


OUT 
OUTPUT 


IN 


“OL HALFD 


Its arguments are: 
XINPUT The variable name of the input quantity. 
HALFD The absolute value of one half of the hysteresis 
band. 


OUTPUT The absolute value of the output of the relay. 


FUNCTION REL3 (XINPUT, OUTPUT) 


FUNCTION REL3 simulates an ideal three position relay. 
OUT 
OUTPUT 


Foch 


IN 











Its arguments are: 
XINPUT The variable name of the input quantity. 


OUTPUT The absolute value of the output of the relay. 


FUNCTION REL3D (XINPUT, HALFD, OUTPUT) 


FUNCTION REL3D simulates a three position relay with dead zone. 
OUT 


OUTPUT 


\ HALFD 


Its arguments are; 
XINPUT The variable name of the input quantity. 
HALFD The absolute value of one half of the dead zone. 


OUTPUT The absolute value of the output of the relay. 


FUNCTION REL3DH1 (XINPUT, DROPOUT, PULLIN, OUTPUT) 
FUNCTION REL3DH1 simulates a three position relay with dead zone 
and hysteresis. Since it involves the storage of past information, it 


can only be used once in a given problem. For this reason, additional 


subprograms, FUNCTION REL3DH2 and FUNCTION REL3DH3 are also available. 
OUT 


OUTPUT 








IN 
“~ PULLIN 


\ DROPOUT 
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Its arguments are: 
X INPUT The variable name of the input quantity. 
DROPOUT The absolute value of input at which the relay 
will drop out. 
PULLIN The absolute value of input at which the relay 
will pull in. 


OUTPUT The absolute value of the output of the relay. 


FUNCTION SATAMP (XINPUT, CUTOFF, GAIN) 


FUNCTION SATAMP simulates a saturating amplifier. 
OUT 





Its arguments are: 
XINPOT The variable name of the input quantity. 
CUTOFF The absolute value of the output when it has 
saturated. 


GAIN The value of the linear gain of the amplifier. 


4, Sample Problem. 

The demonstration problem presented here illustrates the method 
of data presentation required by SIDES. The completed SIDES Coding 
Forms indicate the prescribed method of writing equation statements 
and selecting program options. The actual computer output appears next. 
Both data print-out and single and multiple curve graphs are included in 


this output. 
29 
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